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X-ray reverberation in IH 0707-495 revisited 
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ABSTRACT 

The narrow-line Seyfertl galaxy IH 0707^95 has previously been identified as showing 
time lags between flux variations in the soft- (0.3-1 keV) and medium-energy (1-4 keV) X- 
ray bands that oscillate between positive and negative values as a function of the frequency of 
the mode of variation. Here we measure and analyse the lags also between a harder X-ray band 
(4-7.5 keV) and the soft and medium bands, using existing XMM-Newton data, and demon- 
strate that the entire spectrum of lags, considering both the full energy range, 0.3-7.5 keV, 
and the full frequency range, 10^^ ^ ^ < 10~^ Hz, are inconsistent with previous claims of 
arising as reverberation associated with the inner accretion disk. Instead we demonstrate that 
a simple reverberation model, in which scattering or reflection is present in all X-ray bands, 
explains the full set of lags without requiring any ad hoc explanation for the time lag sign 
changes. The range of time delays required to explain the observed lags extends up to about 
1800 s in the hard band. The results are consistent with reverberation caused by scattering of 
X-rays passing through an absorbing medium whose opacity decreases with increasing en- 
ergy and that partially-covers the source. A high covering factor of absorbing and scattering 
circumnuclear material is inferred. 



Key words: galaxies: active - X-rays: galaxies 
vidual: IH 0707-495 



accretion, accretion disks - galaxies: indi- 



1 INTRODUCTION 

IH 0707-495 (z = 0.0411) is a narrow-line Seyfertl galaxy, 
notable for the presence of an unusually deep dro p in its X- 
ray spectrum just above 7 keV jSoller et al.l l2002h . The spec- 
tral drop has been mode led as an edg e whose energy varies be- 
tween observations ( Gallo et al.l |2004) and the overal l spectral 



form can be explained using partial-cover ing models (B oiler et al 
l2002l:lTanaka et al.l2004l : lGallo et al.l2OO40 . However. iFabian et al 
1 20040 suggested instead that the spectral drop may be the blue 



wing of a rel a tivisti cally-broadened Fe Kq emi ssion line. Re cently, 
iFabian et al.l ( l2009l hereafter F09) and .Zoghbi et al.l ( l20ld . here- 
after ZIO) have claimed that the soft X-ray spectrum of 1H0707- 
495 contains strong, blurred Fe L line emission below 1 keV. F09 
and ZIO explain the flux-dependent spectral variability as due to 
varying contributions from a powerlaw continuum and an ionized 
reflector with high metallicity. To fit the spectrum, the reflected 
emission is extremely relativistically redshifted and blurred, with 
contributions down to 1.23 gravitational radii {GM/(?) from an il- 
luminated disk with reflected emissivity oc (where r is the ra- 
dial coordinate from the black hole). Such a conclusion is surpris- 
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ing, given the extreme black hole spin, the very steep reflectivity 
profile required to obtain detectable emission from so close to the 
black hole, and the implied very small size of illuminating source. 
However, F09 argue that the combination of this spectral modelling 
with an analysis of the source's rapid flux variations does provide 
evidence for emission from a few tens of light-seconds of the event 
horizon of a rapidly-spinning black hole, and hence it is important 
to critically assess that model and test it with further analysis. 

In addition to the spectral modelling, F09 and ZIO searched 
for and detected time lags between flux variations in two X-ray 
energy bands, 0.3-1 keV and 1-4 keV. Lags were analysed in the 
Fourier domain, so that they could be derived as a function of the 
frequency of the mode of flux variation. Those analyses revealed 
the first clear detection of "negative lags" in the X-ray time se- 
ries of an active galactic nucleus (AGN), where negative is de- 
fined as meaning that flux variations in the softer energy band lag 
variations in the harder energy band. Negative lags were found at 
high frequencies (i^ > 6 x 10^* Hz) of amplitude up to about 
40 s. Larger positive lags were found at lower frequencies. Posi- 
tive lags are well-known in AGN (e.g. iMcHardv et al.ll2004h . but 
a pre vious detection of a n egative lag was not strongly signif- 
icant jMcHardvetZll2007h and the measurement of both posi- 
tive and negative lags in an AGN had not been clearly established 
prior to the work of F09 and ZIO. F09/Z10 interpreted the high- 
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Table 1. Summary of observations used: epoch, observation ID, start date, 
duration rounded to tlie nearest ks, and mean 0.3-7.5 keV count rate. 



epocli 


ID 


start date 


duration /ks 


count rate /s ^ 


02 


0148010301 


13-10-2002 


68 


3.841 ± .007 


03 


0506200301 


14-05-2007 


36 


2.028 ± .007 


03 


0506200401 


06-07-2007 


17 


4.227 ± .016 


03 


0506200501 


20-06-2007 


34 


6.116 ± .013 


04 


0511580101 


29-01-2008 


93 


3.504 ± .006 


04 


0511580201 


31-01-2008 


68 


5.163 ± .009 


04 


0511580301 


02-02-2008 


68 


4.509 ± .008 


04 


0511580401 


04-02-2008 


68 


3.656 ± .007 



frequency negative lags qualitatively as a reverberation signal orig- 
inating within a gravitational radius of the black hole event horizon, 
and argued that this interpretation was consistent with their spectral 
model, with its large component of relativistically-blurred reflected 
emission in the soft band. ZIO recognised that the larger-amplitude, 
low-frequency positive lags could not arise from such inner-disk 
reflection, and therefore suggested that these lags might have an 
origin other than reverberation. However, .Miller et al, (2010.) anal- 
ysed the positive lags in the time series of Suzaku X-ray flux varia- 
tions of NGC 405 1 , and recognised that the energy- and frequency- 
dependence of those lags may be well-explained by reverberation 
from substantially more distant material, a few thousand light- 
seconds from the illuminating source, leading to the question of 
whether such a reverberation model is relevant also to 1H0707- 
495. 

In their timing analysis of IH 0707^95, F09 and ZIO did 
not report any measurement of lags for energies above 4keV, 
thereby losing the diagnostic potential of measuring a band that 
contains FeK emission. This is particularly important given that 
the F09/Z10 spectral model predicts that an even larger fraction of 
reflected light should be present in the Fe K spectral region. In this 
paper we show measurements of the lags at higher energy and in- 
vestigate both whether the ZIO model can explain the full set of 
observed lag behaviours and whether a reverberation model such 
as that inferred for NGC 4051 is able to explain the frequency- 
dependent positive and negative lags seen in IH 0707^95. 



2 OBSERVATIONS AND DATA REDUCTION 

XMM-Newton has observed IH 0707^95 during four epochs, 2000 
Oct (Ol), 2002 Oct (02), 2007 Jan-Feb (03) and 2008 Jan (04), 
where the names in parentheses indicate the naming convention 
used by ZIO. As we might worry that different behaviour arises in 
low or high flux states, we only consider observations whose flux is 
comparable to that of the long 04 observations. Thus, in this paper 
we present the detailed analysis of data from 2002-2008, excluding 
the first of the 2007 snapshots (see Table[T]and ZIO). 

The European Photon Imaging Camera (EPIC) observations 
were made with the medium filter and with the prime full win- 
dow mode in 2002 and the prime large window mode in 2007/8. 
Data were processed using SAS vS.O.O and HEAsoft v6.8. Our 
analysis used only the pn data, which offered the best S/N ratio 
of all EPIC data and which were free of pile-up effects. The pn 
events were filtered using standard criteria and instrument patterns 
0-4 were selected. Source events were extracted from circular re- 
gions of radius 35" centred on the source, and background events 
sampled from regions on the same chip about a factor 5 larger in 
area than the source cell. The 2007/8 observations suffered from 
high background, with some background flaring. As noted by ZIO, 



the high background is problematic owing to prominent emission 
lines from the electronics circuit board. The geometry of the cir- 
cuit board results in the cleanest area being the central region of 
the pn chip array, co mprising about the inner 40% of the chip at 
the nominal aim point jLumb et al.l2002h . However, with the target 
positioned in the clean area, it is difficult to find a sufficiently large 
area for similarly clean background extraction on the same pn chip. 
In order to be able to use the largest possible background extraction 
area while minimizing contamination by the background emission 
lines, we applied a tight background screening to the data and re- 
moved periods when the full-band pn background count rate ex- 
ceeded 0.25 count s^^. After screening, the data had a background 
level < 1% of the source count rate in the 0.3-1 keV and 1-4 keV 
bands and 10-15% of the source count rate in the 4-7.5 keV band. 
The total effective pn exposures were 68 ks (02), 87ks(03) and 
297 ks(04). 1H0707-495 ranged in flux from < 1 count s"^ in the 
lowest observed states (not analyzed here) to ~ 10 count s^^ in the 
higher states observed with XMM-Newton. 



3 TIME LAG MEASUREMENT 

F09 and ZIO have previously considered the time lags only between 
two bands: a "soft" band (0.3-1. OkeV) and a "medium" band (1.0- 
4.0 keV). We add a further "hard" band covering 4.0-7.5 keV, where 
the upper energy bound has been chosen to include the Fe K region 
but avoid the region of greatest residual background contamination. 

To m easure time lags b etween these bands we follow the 
method of iMiUer et all ilOldi . which uses a maximum-likelihood 
formalism to measure the best-fitting power spectral densities 
(PSD) in each band and the cross-power spectral densities and time 
lags between bands. To measure a PSD alone, the PSD is defined as 
the Fourier transform of the autocorrelation function, A{t), as in 
other methods, but rather than trying to invert directly the observed 
autocorrelation function, instead we maximise the likelihood, jC, of 
the observed A{t), 



(27r)^/2|C|'/' 



exp 



(1) 



for a time series containing A'' data samples, where C is a A'^ x A'^ 
model covariance matrix, whose elements are given by the model's 
A{t), a is the time series vector of fractional fluctuations with 
length A*' and A"^ is its transpose. The likelihood function as- 
sumes the process that generates the time series to be gaussian 
(or, equivalently, assumes that the Fourier transform of the intrinsic 
time series, before being observationally sampled, has uncorrelated 
phases): in principle this assumption could be varied and replaced 
by a lognormal distribution, however there is no strong evidence in 
this time series for non-gaussian behaviour, and the assumption of 
gaus sianity is common to alternative methods of PSD estimation 
(e.g. lUttlev et aLll2002l) . To find the best-fitting PSD we compute 
the expected covariance matrix for an initial set of model paramete r 
values and iterate using the method of lBond, Jaffe & Knoxl ( Il998l) . 
The model covariance matrix includes the effects of shot noise (see 
iMiller et a l. 2010). The model parameters are the amplitudes of the 
PSD in discrete ranges of Fourier frequency (and in the case of es- 
timation of the cross-spectrum for a pair of time series, the model 
parameters also include the time lags, see below). 

The method has the principal advantages over other ap- 
proaches that it is immune to gaps in the time series data or to 
uneven time sampling with no need for any interpolation (because 
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the likelihood function is only defined for the observed time se- 
ries); it corrects rigorously for the presence of shot noise (includ- 
ing the effect of uncertainty in the shot noise); and it allows rapid 
estimation of the uncertainties in the PSDs and time lags, includ- 
ing estimation of the covariance between frequency bins. Crucially, 
although PSD uncertainties are not gaussian, and neighbouring fre- 
quencies have correlated errors, because the likelihood is calculated 
in the time domain the PSD uncertainties and their covariance are 
straightforward to calculate using either of two methods. First, as 
part of t he likelihood itera tion, we calculate the model Fisher ma- 
trix (see ISond et ^1 19981 for full details) and when the model is 
a good description of the data, the square-root of the diagonal ele- 
ments of the inverse Fisher matrix yield an estimate of the uncer- 
tainties on the PSD parameters. Strictly, these estimates are lower 
limits to the true uncertainties (the Cramer-Rao bound) and do not 
allow for the correlations between parameters which are encoded in 
the off-diagonal elements of the inverse Fisher matrix. To address 
these issues we may instead step through a set of values for each 
parameter and measure the likelihood ratio, £/£niax, with respect 
to the best-fitting model, allowing other parameter values to float in 
the maximisation (a procedure similar to marginalisation over those 
other parameters). In this case, we expect ~ — 21og(£/£max) 
to be distributed as with one degree of freedom. Thus 68 percent 
confidence intervals are found when x^ increases by unity. The 
second method is more time-consuming, of course. The error esti- 
mation mettiodsjiavebeen tested against Monte-Carlo realisations 
for the lMilleretal] j2010l) observation of NGC4051 by Miller (in 
preparation) and found to be in excellent agreement with the distri- 
butions obtained from those realisations. 

To measure ti me delays between two time series, we follow 
iMiller et al.l feOld) and extend the likelihood maximisation to in- 
clude the cross-spectrum of the pair of time series, by including 
both time series in the vector A in equation[T] so that the max- 
imisation process involves handling matrices of rank 2N. The time 
delays are then obt ained from the phases of the cross-spectrum (e.g. 
lNowaketal.ll 19991) . 

Although in principle the multiple observations described in 
Table[T]may be analysed as a single time series, in practice this 
is computationally time-consuming at high time sampling as it 
requires the inversion of several 2A*' x 2N matrices, a process 
which becomes prohibitively expensive on desktop computers for 
2N > 6000. We therefore follow an approach equivalent to that 
of Z 10 in which we ignore the temporal cross-correlation between 
time samples from different orbits, but combine the likelihoods de- 
duced from each observation to obtain an overall best-fit set of PSD 
and time lag values. 

Time series were generated with 48 s sampling. PSDs were 
measured for frequency bins of width Alogj^gZ^ — 0.3 for fre- 
quencies lower than the Nyquist frequency of 0.0104 Hz. The bin 
widths were increased for the lowest frequency bins to ensure that 
their widths were at least twice the lowest observable frequency of 
the longest individual dataset, so that neighbouring frequency bins 
are largely uncorrelated. Fig.[T]shows the PSD estimates in each of 
the three energy bands. The vertical error bars show the expected 
68 percent confidence intervals expected from the square-root of 
the diagonal elements of the inverse Fisher matrix, and hence do 
not allow for covariance between parameters: we do not calculate 
stepped-likelihood uncertainties as we make no further analysis of 
the PSDs in this paper. The PSDs are steep above 10~* Hz, al- 
though the harder bands have flatter PS Ds than the soft band (as 
found also in other AGN, e.g. NGC4051,!MillereLal. 2010). 

Fig.[2] shows the derived time lags as a function of frequency 
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Figure 1. Power spectral densities, uPu, as a function of mode frequency 
u for each of the three X-ray energy bands: soft (open squai'es); medium 
(crosses) and hard (closed circles). Vertical bars indicate 68 percent confi- 
dence uncertainty, horizontal bars show the range of frequencies included 
in each point. Frequency values have been slightly offset between energy 
bands for claiity. 

(the "lag spectrum"). As we now need to sample rapid variations as 
a function of frequency, a higher frequency resolution, A log jq v — 
0.15, was chosen, which also matches more closely the resolution 
adopted by ZIO. Again, the frequency bin widths were increased 
to maintain independence at low frequencies. One further bin at 
lower frequency was included in the maximisation but is not shown 
in Fig.[2]as its uncertainty is too large to be useful. The 68 percent 
confidence region uncertainties shown were calculated by measur- 
ing the likelihood ratio obtained when stepping through parameter 
values, as described above, and thus account for the covariance with 
other parameters and for the non-gaussianity introduced by trans- 
forming from the time domain to the Fourier-phase domain. As we 
have used the combined datasets 02, 03 (high-state) and 04, with 
rejection of high-background segments, and slightly different sam- 
pling from ZIO, we do not expect the medium-soft lags to be iden- 
tical to those of ZIO. There is however good general agreement, 
with maximum positive lags of 237 ± 72 s and maximum negative 
lags of —32 ± 5 s at higher frequencies. The uncertainties appear 
similar to those obtained by ZIO. 

We can see that the hard-soft lags, not reported or discussed 
by F09 or ZIO, follow the same pattern with a larger amplitude of 
time lag. The maximum positive lag in this case is 548 ± 266 s, 
with the most statistically-significant individual positive lag being 
277 ± 62 s. The most negative lag has a value — 55 ± 25 s. The 
existence of several neighbouring lags with similar values further 
increases the statistical significance of both the positive and nega- 
tive lags. The overall significance of the negative lags is discussed 
further below. The lags measured between medium and hard bands 
are not independent of the medium-soft and hard-soft lags and are 
generally in reasonable agreement. There is some disagreement at 
the highest frequencies owing to degeneracies in the fit, highlight- 
ing the difficulty of obtaining a unique measurement at high fre- 
quencies. 



4 DISCUSSION 

4.1 The interpretation of time lag spectra 

Before discussing in some detail the explanation of the observed 
lag spectra, there are some general points to be made. 
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Figure 2. Measured lag spectra (defined to be positive for tlie liarder band lagging the softer band), shown as points with error bars, from the cross-band power 
spectra of the medium v. soft bands (left panel), hard v. soft bands (centre panel) and hard v. medium bands (right panel). Vertical bars indicate 68 percent 
confidence regions from the stepped-likelihood method. Horizontal bars indicate the range of frequencies used for each point. Solid curves show the lag spectra 
expected from top-hat model transfer functions discussed in section |43] 



(i) It has become common-place in X-ray timing analysis to in- 
spect the time lag spectrum and to consider each measured fre- 
quency as being in some sense independent of the others. This ap- 
proach is flawed in the case where lags are caused by reverbera- 
tion, because in general we expect any reverberation signal that is 
restricted to some finite range of time delays to produce lag spec- 
tral signals over the full measurable range of frequencies. The most 
extreme case would be for reverberation seen from a narrow annu- 
lus on an accretion disk viewed face-on: in this case the transfer 
function is a delta-function and its Fourier transform has amplitude 
and phase information distributed over all frequencies. In the case 
of the analysis of IH 0707^95, any claim of a reverberation signal 
at high frequencies should be accompanied by a demonstration of 
the expected time lag signatures at all frequencies. 

(ii) In reverberation models, the measured lag spectrum is a rep- 
resentation of the phases of the Fourier transform of the transfer 
function. Low-frequency modes have time lags in the same sense 
as the reverberation lag (i.e., if the hard band contains more re- 
flected emission than the soft band, the low frequency modes have 
positive lags) but for modes whose time period is less than the re- 
verberation delay, the phase wraps around: phases are evaluated in 
the range —n < (f> ^ n and so any apparent time lag, positive or 
negative, is possible. A corollary statement is that high-frequency 
modes cannot generate large apparent time lags (i.e. any inferred 
time lag at some frequency v must satisfy |r| < 1/2!/) regardless 
of the actual time delay. 

(iii) Because in the X-ray band the signals are likely to be a com- 
bination of reflected, time-delayed emission and directly-viewed 
emission, the measured time delay is always less than the true time 
delay (i.e. the true time delay is d iluted by the addition of the direct 
emission - see iMiller et al.ll2010l) . 

(iv) In cross-correlating two bands and measuring the time de- 
lays between them, there is likely to be reflected emission in 
both bands, and thus the inferred transfer function is the cross- 
correlation of the individual transfer functions of each band. 

(v) The existence of small time lags does not necessarily indi- 



cate a short light travel-time between the source and the reflector, 
because reflection from material along the line of sight, at any dis- 
tance from the source, has zero time delay with respect to directly- 
viewed emission, and slightly off-axis material can result in arbi- 
trarily small delays. As an example, consider the transfer function 
of a thin spherical shell of reflection, of radius r, isotropically sur- 
rounding a source. The transfer funct ion is a uniform top-hat distri- 
bution with time lags ^ r ^ 2r/c ( |Petersonlll993h : the existence 
of lags close to zero do not imply that the reflecting material has to 
be close to the source, and such a conclusion could not be reached 
unless a full model reproducing the time lags for all modal frequen- 
cies were to be tested against the data. 

Bearing these points in mind, we shall first critically review 
the previous analyses of the lag spectra of IH 0707-495, and then 
discuss a more viable model that explains all the observed features. 



4.2 The inconsistency with inner-disk reflection 

The observed time lags are inconsistent with the hypothesis that 
they arise from reflection from the inner accretion disk, as pro- 
posed by F09 and ZIO, for a number of reasons. In this section we 
first discuss the hypothesis of F09/Z10, that there are two indepen- 
dent lag-generating mechanisms that operate at different frequen- 
cies, and we then move on to test consistency of the new hard-band 
lag measurements with the F09/Z10 model. 



4.2.1 Constraints on two mechanisms of lag generation 

Consider first the time lags discussed by F09 and ZIO, between 
the soft and medium bands. At low frequencies, v < 10~^'^ Hz, 
the medium band lags the soft band, which is inconsistent with 
the soft band containing delayed reflection of the medium band. 
The signature of inner-disk reflection is supposed instead to be the 
small (~ 30 s) negative lags at high frequencies (v > 10~^ Hz). 
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This was claimed to place the reflecting region a few tens of light- 
seconds away from the Fe L ionising source, traced by the medium 
band. But if there are longer-timescale fluctuations in that source, 
as observed, these must also be delayed in the same sense: i.e. there 
should be a negative lag of tens of seconds at low frequencies. This 
is not observed, instead positive lags of ~ 200 s are seen. F09 
and ZIO suggested these arise as a separate process, such as the 
propagation of fluctuations through an accretion disk from outer 
to inner radii, with the inner disk emission supposedly having a 
harder powerlaw spectrum than the oute r disk, so that the softe r 
band varies sooner than the harder band dArevalo & Uttlevll2006h . 
This scenario could only produce the observed lag spectra if (i) in- 
ward radial perturbations leading to variable emission existed and 
(ii) radial dependence of the spectrum existed and (iii) the softer 
outer disk regions did not vary on short timescales, whereas the 
harder inner disk regions did, and (iv) the outer disk emission did 
not result in any soft-band reflection, implying the outer disk spec- 
trum to be strongly cut off just at the Fe L ionisation edge, and (v) 
despite having emission from a wide range of radii, nonetheless the 
reflected emission had a radial dependence (F07, ZIO). Such a 
set of requirements is not necessary if a more straightforward ex- 
planation is adopted (section l431 >. These considerations motivate us 
to carry out further evaluation of the viability of the ZIO model. 

4.2.2 Testing the model for consistency with the hard-band lags 

The hard band (4-7.5 keV) analysed here also shows significant 
time lags with respect to the soft band, again in the sense that the 
hard band lags the soft band by up to ~ 550 s at low frequencies, 
with negative lags at high frequencies. However, the hard band 
contains the FeK region and, in the Z10/F09 model, should also 
show delays with respect to its ionising continuum. According to 
the model, most reflection occurs in the hard band, then the soft 
band, then the medium band with the least (but non-zero) reflec- 
tion (see Fig. 8 of ZIO). We find that the hard-soft lag spectrum has 
an almost identical form to that of the medium-soft lag spectrum, 
but with lags of larger amplitude, despite the expectation from the 
ZIO model that the hard band supposedly has more reflection than 
the soft band and thus should lag the soft and medium bands at 
v > 10"'^ Hz (i.e. the hard-soft and hard-medium lags should be 
positive where the medium-soft lags are negative). 

We have assessed the statistical significance of the disagree- 
ment with the model as follows. Starting with the medium-soft lag 
spectrum, we find the mean lag of the six frequency bins that have 
negative lags, i.e. 0.00066 < u < 0.0053 Hz, by fitting a a single 
time lag value to that frequency range and evaluating its uncertainty 
using the maximum-likelihood method. All other bandpower am- 
plitudes and time lags were allowed to vary during the likelihood 
maximisation. The 68 percent confidence interval for that one pa- 
rameter is given by Ax^ ~ — 2 log(£/£max) = 1, which yields 
a mean lag rmcdium-soft = —16.4 ± 2.4 s. The mean lags for the 
same frequency range in the other comparisons between wavebands 
are Thard-soft = -37.0±9.0sandrhard-modmm = -14.7±8.7s 
(note that these three lag values are not statistically independent but 
are consistent with each other). Thus all three lag spectra are sig- 
nificantly negative in the frequency range of interest, counter to the 
expectation of the ZIO model. 

To compare with the ZIO spectral model, we estimate the ex- 
pected relation between the lags in the three bands as follows. In 
a waveband where both direct and reflected light is present, with 
reflection fraction /, the inferred lag, r', is diluted with respect 
to the intrinsic time lag to the reverberating region, r, such that 



r' ~ /t/(1 + /) for gir -C tt when considering Fourier modes 
with angular frequency oj jMiller etal. 2010). Thus the inferred lag 
between two such bands, designated i and j, with reverberation sig- 
nals from the same region, is given by 

Ar/, ~ ( r. (2) 

We obtained values for the reflection fractions /soft, /medium, 
/hard in each of the soft, medium and hard bands respectively 
by refitting to the data the model described by ZIO, with re- 
flecti on spectrum given b y relativistically-blurred REFLIONX ta- 
bles URoss & Fabianll2005ln . As the time series analysis was car- 
ried out on variations in count rate and not flux (as this opti- 
mises the signal-to-noise in the timing analysis) the fractions / 
have been obtained as ratios of counts in the blurred reflection 
component to the counts in the power-law component in each 
band. The reflection fractions in the summed observations were 
/soft = 1.60, /medium = 0.57 and /i^ard = 2.03 in the soft, 
medium and hard bands respectively. The formal statistical un- 
certainties on the model fits are very small, however we expect 
there to be some possible variation in allowable models, and to al- 
low for this modelling uncertainty we also fitted the ZIO model 
to the data from each of the four separate orbits that make up 
the main 04 observation. The sets of four reflection fractions 
were found to be /soft = {1.63,1.13,1.57,2.48}, /medium = 
{0.58, 0.40, 0.57, 0.91}, /hard = {2.04, 1.45, 2.05, 3.58} in each 
band for each of the observation IDs 0511580I0I, 0511580201, 
0511580301, 0511580401 respectively. Inserting these reflection 
fractions / into equation|2l we can predict the expected relation 
between the observed lags in each band for any assumed value of 
distance of the reflector from the primary source. 

Fig.[3] shows these expected relations and contours of Ax^ 
obtained from the maximum-likelihood fitting, comparing in turn 
the hard-soft lags and the hard-medium lags with the medium- 
soft lags. Confidence region contours are shown at Ax^ = 
2.3, 4.6, 9.2, 13.8 which correspond to the 68, 90, 99, 99.9 percent 
confidence regions for two parameters. The failure of the hard-soft 
lags to have the opposite sign to the medium-soft lags rules out the 
ZIO model at > 99.9 percent confidence. The medium-soft lags 
also rule out the model at > 99 percent confidence, although these 
two measures are not independent. 

4.3 The reverberation transfer function 

We can use the considerations of section|4T|together with the new 
information from the hard band to infer what reverberation model 
could explain the full set of available information. Given that a re- 
verberation explanation should also seek to explain the large posi- 
tive lags measured at low frequencies, we should identify the form 
of transfer function(s) that would produce the overall lag spectra. 

It is important to note that the lag spectra derive from the 
phases of the Fourier transform of the transfer function. The am- 
plitude information is carried in the cross-band power spectrum. 



^ We note that the model is not a good fit to the data over the entire energy 
range, in particular there are significant residuals of 20 percent in the Fe K 
emission region, as seen in Fig. 8 of ZIO, in addition to significant residuals 
below 5keV that ZIO argue might be an effect of additional variations in 
ionisation. The extreme blurring in the model has erased features ai'ound 
FeK that are present in the data. However, the ZIO model requires such 
extreme blurring, so we retain that model in our test. 
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Figure 3. Confidence regions for the mean lags in the frequency range 
0.00066 < u < 0.0053 Hz between (left) hard-soft and medium-soft 
bands and (right) hard-medium and medium-soft bands. Contours are 
shown corresponding to the 68%. 90%, 99% and 99.9% confidence regions. 
Straight lines indicate the relationship between the lags in these bands ex- 
pected from the ZIO model for each of the four 04 observations (only neg- 
ative medium-soft lags are allowed in the ZIO model). 

but it is mixed in with the unknown underlying cross-band spec- 
trum, and so we have to rely on the phase information alone to 
inform us about the transfer function. Thus we cannot uniquely cre- 
ate a transfer function from the lag information. We can however 
consider some simple models that may have some link to physical 
reality and that create the observed lag spectrum features. 

There are three effects that could lead to lag spectra of the 
form seen in Fig.[2] First, transfer functions with sharp transitions 
in general have oscillatory phases, and if the transfer function has a 
maximum away from zero lag, the phase oscillations can go nega- 
tive. Such transfer functions may be generated by clumpy reverber- 
ating material or m aterial not isotropically surrounding the source 
( lMilleretalJl2010h . As an example, consider a transfer function 
that is a top hat in time delay, of width At centred on delay io and 
with intensity some fraction / of the directly-viewed intensity. The 
lag spectrum, t{uj) for angular frequency oj, is given by 

fsinc(a;At/2) sincjin 

tancjT = \ / , — (3) 

1 + /sinc(ajAt/2) coscjto 

An example is shown in Fig.|4^ for /S.t = to = 1000 s, / = 1. Neg- 
ative time lags are seen at high frequencies, although it is difficult 
to obtain negative lags over a wide frequency range from reverber- 
ation in one band only. However, in section l^Tl fiv) we noted that 
reverberation is likely to be present at differing levels in both bands 
being cross-correlated. Fig.|4j5 shows the lag spectrum expected for 
transfer functions that are uniform from delay t — to ti where 
ti = 1000 and 200 s in each of two bands, and with reflection 
fractions f — 1 and 0.5. Now a broader range of negative frequen- 
cies is seen. Finally we consider the effect of the transfer functions 
having no reflection less than some cutoff: Fig.|4j; shows the same 
transfer functions as (b) but with no reflection at delays t < 100 s. 
This allows additional positive oscillations at high frequencies. 

We can approximately fit this simple top-hat model to the mea- 
surements (Fig.[2j. The models shown were fitted jointly to the 
three lag spectra in two stages. For speed, an initial fitting pro- 
cedure minimised ignoring the covariance between parameters 
and the non-independence of the lag spectra. Because the lag uncer- 
tainties are non-gaussian, and because there is covariance between 
both lag values and cross-band power spectrum amplitudes, to es- 




10-* 10-" 0.01 

v/Hz 



Figure 4. Models illustrating oscillatory lag spectra: (a) thin solid line top- 
hat transfer function in one band only: (b) thick solid line top-hat transfer 
functions in two bands, both transfer functions continuous from zero lag; 
(c) dotted line as (b) but with no reflection delays shorter than 100 s (see 
text for details). 

timate the goodness-of-fit we then took this model and refitted the 
time series using the maximum likelihood method, but with time 
lags fixed to the model values. From this, the goodness-of-fit was 
evaluated as = —2 log(£//Imax) where C was the model like- 
lihood and £max was the maximum likelihood found by allowing 
each frequency bin to have an independent time lag (i.e. the values 
shown as discrete values in Fig.|2ll. Considering just the hard-soft 
and medium-soft lags, to avoid the problem that all three lag spec- 
tra are not statistically independent, the model shown had ~ 23 
for 25 degrees of freedom, a statistically-acceptable fit. The model 
fitted is one of the simplest possible models that we could have 
chosen, and the key point of this exercise is that even this simple 
model generates well the large positive lags at low frequencies, the 
sharp transition, and the extended range of high frequencies over 
which the lags appear small but negative. No such model has been 
demonstrated for the inner-disk reflection hypothesis of F09/Z10, 
which was ruled out by the considerations of section l4".2.2l 

The best-fit top-hat start time was 100 s, assumed common to 
all three bands. The maximum time delays for the hard, medium 
and soft bands were 1800, 1450 and 150 s respectively, and the ra- 
tios / of flux in the reflected component to flux in the directly- 
viewed component (which may itself be partially obscured) were 
0.67, 0.38 and 0.30 respectively. We do not evaluate uncertainties 
on these parameter values as this is computationally extremely ex- 
pensive, requiring stepping through all parameter values and re- 
maximising the likelihood in the time domain, and the model is 
almost certainly too simple to be taken too literally. However, the 
values do indicate a possible physical origin for the different trans- 
fer functions in the three bands. The maximum time delay increases 
with energy: this is expected if X-rays are scattered while travel- 
ling through a medium whose opacity decreases with increasing 
energy, because then in the soft band X-ray photons cannot travel 
far through the medium before being absorbed, and any scattered 
X-rays must have relatively short time delays, whereas hard X- 
ray photons can travel further and create reverberation signals with 
longer delays. In this case we would expect also to see a greater 
fraction of scattered X-rays in the hard band, qualitatively in agree- 
ment with the inferred fractions. However, the observed fractions 
are also affected by any absorption of the directly-seen radiation, 
which cannot be determined from this analysis, and the predicted 
fraction of scattered light depends on the unknown radial distribu- 
tion of scattering and absorbing material, making a more quantita- 
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tive test difficult to achieve. This inference imphes that the rever- 
beration transfer functions appear to arise from reflection or scatter- 
ing from a distribution of material rather than from an unobscured 
reflector such as a naked accretion disk. Given the high fractions of 
reflected light, we infer that either the reverberating material sub- 
tends a large solid angle at the source, or that the directly-viewed 
emission is partially-obscured, or both. 

According to this simple model, the reverberation in the soft 
band is characterised primarily by small delays, 100 < t < 150 s. 
However, as noted in section |4n the existence of small lags cannot 
be taken as evidence for a short light travel time between source 
and reflector. In the case of IH 0707-495, the model is only vi- 
able if the medium and hard bands have time delays up to approxi- 
mately 1800 s. Thus, even if one wished to assert that the soft-band 
reverberation were caused by inner-disk reflection, one would re- 
quire the medium and hard band reverberation to be produced at 
significantly larger distances, with no physical explanation for that 
difference. Such an explanation would not allow models with these 
transfer functions in which both soft-band Fe L and hard-band Fe K 
emission were reflected from the same inner disk region. 

4.4 Spectral modelling and further work 

The aim of this letter has been to investigate the timing information 
for IH 0707-495 rather than to address the spectral model. The 
above analysis strongly indicates that there is a substantial com- 
ponent of scattered light present in the spectrum. Spectral mod- 
elling of this source to date has only considered either extreme 
relativistically-blurred models with little obscuration (F09, ZIO) 
or partially-covering absorption models l lGallo etalj|200 4l) with 
no scattering. ZIO claim to have ruled out the partial covering 
models based on the weak observed strength of Fe Ka emission, 
but we note that t heir argument has already been discounted by 
iMiller et al.l ( l2009h and lYagoob et alj j201(t) . In fact, both sets of 
published models require unattractively high values of Fe abun- 
dance, around 1 1 for the relativisticall y-blurred model ( F09) and 
around 3 for the partial-covering model dOallo et al.ll2004l) . The in- 
dication that there is a substantial fraction of X-rays scattered on 
passing through an absorbing medium means that models that in- 
clude both absorption and scattering are required. The smooth spec- 
tral features in the soft band indicate that the scattered emission is 
likely to be smoothed by Compton scattering as seen in the Monte- 
C arlo radiative tran s fer mo dels of transmission through a disk wind 
of ISim et al. U2008ll20I0al lbl). and the presence of a large amount 
of scattered light provides a natural way of obtaining a high equiv- 
alent width in a broad, smooth emission feature, which is likely 
combi ned with st rong absorption edges. Further modelling such as 
that of lSim et"^ (oD. cit.) is required to make progress. 



5 CONCLUSIONS 

We have measured and analysed the time lags in I H 0707^95 
between three X-ray bands: the soft (0.3-1 keV) and medium (I- 
4keV) bands discussed by F09 and ZIO, and the hard band (4- 
7.5 keV) not previously considered by those authors. Time lags 
have been analysed in the Fourier domain. We can summarise our 
conclusions as follows. 

(i) The F09/ZI0 model requires separate mechanisms for the 
positive lags at low frequencies and the negative lags at high fre- 
quencies, and in Section l4.2.II we have argued that this imposes 



some special constraints on those mechanisms which may not be 
fulfilled in practice. 

(ii) Consideration of the lags in all bands, soft, medium and 
hard, rules out models in which the inner accretion disk makes a 
significant contribution at both FeL and FeK energies, and thus 
rules out the spectral model of F09/Z10. In particular, the lags be- 
tween the hard and medium bands and between the hard and soft 
bands are also negative in the same range of frequencies that the 
lags between medium and soft band are negative. Comparison be- 
tween the observed lags and the lags predicted from the spectral 
model rules out the ZIO model at confidence > 99.9 percent. 

(iii) A simple model of reverberation, in which all X-ray bands 
have differing levels of reverberation, provides a good quantitative 
description of the large positive lags at low frequencies and the high 
frequency transition to negative lags, and fits the full lag spectra 
between all bands given the measurement uncertainties. 

(iv) The harder X-ray bands appear to require longer maximum 
time lags than softer bands, consistent with a geometry in which re- 
verberation arises, not as reflection from a naked accretion disk, but 
from scattering of radiation that is passing through partially-opaque 
material. Time delays at least as large as 1800 s are required, plac- 
ing a lower limit on the extent of the reverberating region of approx- 
imately 1000 light-seconds, corresponding to 20 -100 gravitatio nal 
radii for black hol e mass i n the range 10^ M© jLeighlvll2004l) to 
2 X lO'^ Mq (ZlO. IZhou& Wang 2005) respectively. 

(v) Spectral models that include the effects of both absorption 
and Compton scattering are required. 

These results strengthen previous conclusions that X-ray observa- 
tions of the inner regions of type-I AGN may be dominated by the 
effects of circumnuclear material a few light-hours from t he cen- 
tral so urce (e.g. Miller e t al. 2008; Risaliti et al. 2009; Mil ler et all 
I2OIOI) . The material causes both scattering and absorption, and 
subtends a large solid angle at the source but with an anisotropic 
geometry. Coupled with the evidence for e nergetic outflows (e.g. 
[Pounds & R eeves''2009l; lReeves et ai]|2009l) it is likely that the re- 
verberating material is associated with an accretion disk wind. 
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